Higher dimensional black holes with a generalized gravitational action 

Jerzy Matyjasel^ and Malgorzata Telecka 
Institute of Physics, Maria Curie- Sklodowska University 
pi. Marii Curie- Sklodowskiej 1, 20-031 Lublin, Poland 

Dariusz Tryniecki 

Institute of Theoretical Physics, Wroclaw University 
pi. M. Borna 9, 50-204 Wroclaw, Poland 
(Dated: February 1, 2008) 

We consider the most general higher order corrections to the pure gravity action in D 
dimensions constructed from the basis of the curvature monomial invariants of order 4 and 
6, and degree 2 and 3, respectively. Perturbatively solving the resulting sixth-order equations 
we analyze the influence of the corrections upon a static and spherically symmetric back hole. 
Treating the total mass of the system as the boundary condition we calculate location of the 
event horizon, modifications to its temperature and the entropy. The entropy is calculated 
by integrating the local geometric term constructed from the derivative of the Lagrangian 
with respect to the Riemann tensor over a spacelike section of the event horizon. It is 
demonstrated that identical result can be obtained by integration of the first law of the 
black hole thermodynamics with a suitable choice of the integration constant. We show that 
reducing coefficients to the Lovelock combination, the approximate expression describing 
entropy becomes exact. Finally, we briefly discuss the problem of field redefinition and 
analyze consequences of a different choice of the boundary conditions in which the integration 
constant is related to the exact location of the event horizon and thus to the horizon defined 
mass. 

PACS numbers: 04.50.+h, 04.70.Dy 



I. INTRODUCTION 

In recent years gravitation theories with higher derivative terms have attracted a great deal of 
attention. Indeed, according to our present understanding the general relativity is to be treated 
as the lowest order term of the effective theory consisting of a series of loop or string corrections. 
Typically such corrections are constructed from higher powers of curvature and their derivatives, 
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and, hence, the gravitational action I can be written as 

m 

I = '^aklk, (1) 

k=0 

where for > 1 is the contribution of operators of dimension 2k, Iq is related to the cosmo- 
lo^cal constant and ..e arbitrary constant. Anron^the hrgher cnrvatute theories a great deal 
of activity has been focused on the Lovelock gravity 1]. In this theory, the Lag rangian Lm is 
constructed from the dimensionally extended Euler densities of a 2A;-dimensional manifold 

r n-fc i:Ciiii...Cfedfe p aibi p aj.fefc /r,\ 

J-k — ^ "aib-i_...akbk^cidi ■■■^Ckdk ' V^/ 

where the generalized 5 function is totally antisymmetric in both sets of indices. A m— th order 
Lovelock action, 1^, is the sum of m + 1 terms given by Eq. ((2j of ascending complexity 

Im = I d''x{-g)^I^L^ = j d^x{-gfl^Y.auCk, (3) 

fc=0 

where at are arbitrary parameters. If the Lovelock action includes the terms up to Cm, the 
dimension of the spacetime should satisfy D > 2m + 1. 

Varying the action functional I^a with respect to the metric tensor one obtains the equations of 
motion of pure gravity, which can be schematically written as 

Originally, Lab has been introduced by Lovelock to demonstrate the most general symmetric and 
divergence-free second rank tensor, which can be constructed from the metric and its first and 
second derivatives. Since the resulting equations of the Lovelock gravity involve at most second 

?her curvature theories 
which do not 

approach their Einsteinian counterparts as the perturbative expansion parameter is set to zero, 
and, when linearized, the Lovelock equations are free of ghosts. Moreover, the higher-order terms 
appear quite naturally as the low- energy limit of the string theory P, . 

At each order Ck is a linear combination of the basis curvature invariants with the particular 
set of coefhcients calculated from Eq. ©. For example, in the first two (simplest) cases, one has 
a cosmological constant and the curvature scalar, for A; = and k = 1 respectively. At A; = 2 there 
are three invariants which are combined into the Gauss-Bonnet term whereas at k = 3 the basis 
has eight members. 



derivatives of the metric it avoids some of the typical problems of the hig 
P, P . Specifically, at the classical level, it avoids singular perturbations J. 
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On the other hand, one may consider the more general curvature terms, with arbitrary coeffi 



cients rather than those inspired by the particular theory. (See for example [3, |a| and references 
cited therein). In this case, the relation between dimension of the spacetime and the order of higher 
curvature terms retained in the action functional is lost and the dynamical equations inevitably 
involve higher derivatives of the metric. There is nothing wrong in using such equations, pro- 
vided only the physical solutions are selected. However, the expected complexity of the resulting 
equations may be a serious obstacle in this regard. 

The literature devoted to various aspects of the higher derivative gravity is vast indeed. As the 
examples of such theories one may consider the quadratic or higher order gravity (see 

and the references cited therein), and, when a negative cosmological constant 
is present, the Einstein-Hilbert action in D = 5, supplemented with Riem^ = Rabcd.R'^'"^'^ : which 
corresponds to an effective AdSs bulk action Ql- 

The foregoing discussion indicates that the analyses can be carried out in two directions. First, 
one can construct and investigate the possible candidate terms that may appear in J^, whereas the 
second direction of calculations is to apply the thus constructed equations in the particular context 
of black hole physics or cosmology with or without additional matter fields. The latt er appr oach 
has been successfully applied in various contexts in Refs. 
in the references cited therein. 
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A few words of comment are in order here. First, it should be observed that we have no 
information on m, i. e., the number of terms that should be retained in Eq.(^. Second, and that is 
closely related to the above observation, it is really desirable and perhaps unavoidable to construct 
the solutions of the dynamical equations which reflect the nature of their derivation. Finally, 
it should be observed that because of complexity of the problem the full system of equations is 
probably intractable analytically and one has to construct either approximate solutions or refer to 
the numerical methods. 

In this note we shall explore the second possibility and perturbatively solve the equations 
resulting from the variation of the D— dimensional action without a cosmological constant 



about the Tangherlini black hole 



271 . l28l |. In doing so the class of solutions are restricted to 



the admissible ones. We shall assume that the total action functional / is the sum of the flrst 
three nonvanishing terms, 1^, constructed from the basis curvature invariants. That is, we assume 
arbitrary coefficients rather than those inspired by the particular theory. The results of this paper 



generalize results of Lu and Wise 



and may be though of as a partial genera 



analogous results obtained within the framework of the Lovelock gravity 



;ial generanzation 
IE, 30, 31, 3- 



ization of the 
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II. EQUATIONS 



We shall consider the action functional being a sum of the terms (conventions are Rab = R'^acb ~ 



dcTl;^ , signature -,+,+,+) 



h = al d''x{-g)^/^R, (5) 
l2 = J d^x{-gY'^ (h^R^ + b^RabR''' + b^R^.^aR'"'"') (6) 



and 



h = j d'^xi-g)^'^ (ciR" + C2RRlRl + csRRj'^Rj' + c^RlRtRS" 

+ C^R^RbR-c + (^^R^Rbc'' Rde" + C-lRab'^Rj^ ^ef"^ + ^sRce'^^Raf'^Rbd'^) ' (''') 

where a = (IGttGd)"^ and G_d is Newton's constant. That is, we will restrict ourselves to scalar 

n 

terms of order 2, 4 and 6 belonging to classes i! ^4 2 ^6 3' respectively ^Sj. In the course of 
the calculations we shall assume that the coefficients bi and Ck satisfy \bi \ / a « 1 and \ck/bi\ « 1 
for i = 1, .., 3 and k = 1, 8, respectively. The case bi ^ Ci can be easily obtained from the former 
one simply by relegating the terms proportional to bibj from the resulting expressions. 

It should be noted that depending on the dimension D the curvature terms may be subjected to 
additional relations Q]. Moreover, for a static and spherically symmetric line element we have addi- 
tional vanishing combination of the elements of the curvature basis with the coefficients depending 
on D 



L vaij 



Although it is possible to adopt (with small modificatioiis) the results presented in Refs. 



l35l | , here we proceed differently and use the Weyl method 



36, 



37 



33, 



34 



38| | . The line element describing 



the static, spherically symmetric D = d + 2-dimensional geometry may by cast into the form 

ds^ = -f{r)dt^ + h~^ir)dr^ + ^ ''f^ r^dx'dx^ i,j = 2,...,d+l, (8) 

1 — x"^ 

where Xi are the coordinates covering maximally symmetric d -dimensional space. The components 
of Riemann tensor, the basic ingredients of our calculations, are simply 

i?,/- = r-\l - h') (^f 5f - 6r5^) , (9) 



i?^/'^ = -r-i/i/i'<5^', (10) 
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Rtr' = -f-'hifh)' (11) 

and 

Ru'' = -r-'fih^f'Si, (12) 

where the prime denotes differentiation with respect to the coordinate r. Upon inserting the line 
element into / and subsequently varying the thus obtained reduced action with respect to the 
functions / and h, one obtains rather complicated system of equations (not displayed here) , which 
may be further simplified with the substitution 

= (13) 

and 

Except for certain combinations of the numerical coefficients bi and Cfe leading to the dimensionally 
extended Euler densities, the resulting equations of motion are still too complicated to be solved 
exactly. Fortunately, one can easily devise the perturbative approach to the problem, treating the 
higher derivative terms as small perturbations. 

Now, in order to simplify calculations and to keep control of the order of terms in complicated 
series expansions, we shall introduce another (dimensionless) parameter e, substituting bi — > ebi 
and Ci s^Ci- We shall put e = 1 in the final stage of calculations. For the unknown functions 
M(r) and ip{r) we assume that they can be expanded as 

m 

M{r) = e'Mi{r) + 0(£"*+^) (15) 

i=0 

and 

m 

^P{r) = Ye'Mr) + Oie^^'). (16) 

i=l 

The system of differential equations Mi{r) and ipi{r) is to be supplemented with the appropriate, 
physically motivated boundary conditions. First, it seems natural to demand 

M(r+) = (17) 
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or, equivalently, Mq (r+) = and Mj (r+) = for z > 1, where r+ denotes the exact location 

of the event horizon. Such a choice leads naturally to the horizon defined mass, 

On the other hand, however, one can relate the additive integration constant with the total mass 
of the system as seen by a distant observer. Indeed, analysis carried out in a weak field regime 
indicates that the constant of integration C is related to the mass M according to the formula 

(^^) 

where 

27r(d+i)/2 

"'^ = m+i)/2) ^^'^ 

is the area of a unit d-sphere. For the function tp{r) we shall always adopt the natural condition 
^(oo) = 0. Since the results obtained for each set of boundary conditions are not independent and 
one can easily transform solution of the first type into the solution of the second type (and vice 
versa), we shall concentrate on the boundary conditions of the second type. A brief discussion of 
the consequences of the boundary conditions of the first type will be given at the end of the paper. 

III. PERTURB ATI VE SOLUTION 

One expects, on general grounds, that the terms proportional to the coefficients ci, C2, C4 and 
C5 do not contribute to the solution. Integration of the zeroth-order equations yields, as expected, 
the Tangherlini solution with Mo{r) = C, whereas the first-order equations give 



^^^(r) = -^{d-2)id-l) (21) 



and V'l = 0. Integration of the second-order equations, although straightforward, yields much more 
complicated results: 

-^2(r) = ^l^^2+2d'' [4c3ci(5 + 8d + 3^2) - 2c6(3 - 5d^ - 2d^) - 2c7(10 + 16d - 21^^ - 7d^) 

+ ^C8(2 + 2bd - 24d^ - 7d^) - ^^{d -2){d + l){3d + 5)- ^^^(d -2){d + l){d + 2) 



oz,2 

■^(d-2)(4d2 + 19d + 9) 
a 



C\d'-1) jd + l) 

(lj-d+3 



_ 3c»«i - 1, - (rf _ 2) _ 1?«L(, _ 2) _ _ 2) 



[Scad - C6(2 - 3d) + 12c7(d - 1) 

(22) 
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and 

'^2(^) = ^^2d+2 [4c3d(3 + 2d) - 2c6(l -2d- d^) - 6c7(2 -2d- d^) 

+^C8(2 - 3(i - d^) _ ^-^{d - 2){2d + 3) - -2){d + 2)-^{d- 2){2d + 5) 

2 a a a 

(23) 

If 6j ~ Ck then the terms involving hi in M2(r) and V'2(^") should be omitted. Note that for Z) = 4 
{d = 2) the function Mi(r) is identically zero, whereas in M2{r) and 'ip2{i') all the terms proportional 
to bihj = 1, 2, 3) are absent. It is because the Kretschmann scalar, Rabcd.R"''"^'^ , can be relegated 
from the action as the functional derivative of the Gauss-Bonnet term with respect to the metric 
is zero. 

The line element ((HJ with H2H23() provides the most general solution to the problem. Since 
the partial results referring to a particular dimension and/or the concrete form of the action exist 
in the literature, it is worthwhile to compare them with our solution. For D = 4 our solution 
reduces to that which can be easily constructed by integrating equations derived by Lu and Wise 
[2^. Similarly, retaining in the action functional Q only the term proportional to cy and making 
substitutions C = M/Ml, a = MI/IQ-k and C7 = a/Mp, where Mp is the Planck mass and a is 

r 1 

a coupling constant, one obtains results presented in Ref. 21], however, with one reservation: We 
have not observed any difference between the values of the radii at which qqq and g"^^ vanish, that 
is, of course, in concord with the general form of the adopted line element. 

Before proceeding further let us compare our solutions to the analogous solutions of the Lovelock 
gravity. First, let us introduce a = ai, bi = a2^i and Cj = a^Ci, where Q2 and Q3 are the numerical 
coefficients that stand in front of the Gauss-Bonnet and the third-order term of the Lovelock 
Lagrangian, whereas 6j and are the coefficients that give rise to the Lovelock terms. Numerically 
one has 61 = 63 = —62/4 = 1 and ci = 1, 62 = — 12, C3 = 3, C4 = 24, C5 = 16, cg = —24, cy = 2, cg = 
—8. Inserting the above coefficients into the solution, after massive simplifications, one obtains 
fi^ (r) = (r) (or, equivalently ^(r) = 0) and 



h\r) = 1 - ^ + 4096 02 



l_512^?^|^(d-2)(d-l) 



{d-2){d-l) 



7^4(74 A^3 

+ «365536 ^^d+M - 4)(d - 3)(d - 2){d - 1). (24) 
On the other hand, although the action of the Lovelock gravity looks rather complicated it is 



29 



possible to construct an exact solution describing static and spherically symmetric configuration 
30, 3ll- 

Such a solution can be found, after the substitution f'^ = h? = l — r'^F{r), by solving for 
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real roots of the m — th order polynomial equation 



m 



E^^F' = ^ (25) 

k=0 ' 



where the coefficients q are given by 



"0 1 - rn(i\ 



and 

2k 



Cfc = — JJ((i + 2-n) for A; > 1. (27) 



n=6 



Note that if the cosmological constant is taken to be zero then qq = 0. Finally, assuming a2 = 6 ~ e 
and 03 = c ~ and expanding Eq. (|25() in powers of e, one easily reproduces Eq. ()24() . 



IV. TEMPERATURE AND ENTROPY 



The approximate location of the event horizon of the black hole solution derived in the previous 
section is given by 

r+ = (2C)i/('^-i) - ^63(2C)-i/('^-i) - -(2C7)-3/('^-i)|c3(d2 _ i)d + -ce{d^ - l){d - 1) 
a a 2 

- Ic7[2 -{d- l)(5d - 4)d] + ^C8(10 - 13d + 12(i2 - 5d^) - -(6163 + &263)(d^ - '^){d - 2) 
2 8 a 

+ ^bl{d-2){4 + 6d-13d)}. (28) 
For D = 4, the first order terms are absent and one easily reproduces Lu-Wise result 



2GM 



(29) 



As is well-known, the Hawking temperature, T, can be easily calculated from the metric tensor 
without referring to the field equations. The standard method of obtaining T relies on the Wick 
rotation. The Euclidean line element has no conical singularity provided the time coordinate is 
periodic with a period P given by 

P = 47r ^hm {gttgrrf^ • (30) 

Its reciprocal is identified with the black hole temperature, which, in the case on hand, reads 

T = ^(2C)-V{^-i) _ (rf-2)(d-l)rf^ 3/(.-i) 
Air A-Ka 

+ l^(2^)"'^^'"'^i-|('^ - 4)[2 + (d - 2){d - l)d] 

+ ^{d - 4)[4 + (d - 5)(d - l)d\ - ${d - 2)2(4 + 7d- 4d^)}. (31) 
o 2(1 
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It should be noted that the Hawking temperature does not depend on 61, b2, C3 and cq. This 
behavior can be traced back to the possibihty of removing curvature terms proportional to this 
very coefficients by means of the appropriate field redefinition. Such a redefinition certainly modifies 
equations of motion of the test particles but should not modify the temperature, which, in turn, is 
to be measured at infinity. On the other hand, the horizon defined mass leads to the expression for 
temperature which depends on the full set of parameters. The reason is that the horizon defined 
mass is not the mass measured at infinity. We shall return to this problem later. 

From Eq. pij) one sees that the heat capacity Cbh = OM/dT calculated for the Tangherlini 
black hole is given by 

and is always negative. That means that for d > 1 such a black hole is thermodynamically 
unstable, i.e., it increases its temperature when radiating. On the other hand, for T >> 1 the 
higher derivative corrections can modify this behaviour, and, depending on the signs and values of 
the coupling constants they can give a nonnegative contribution to the total heat capacity as can 
be easily inferred from 

Cbh = C]^H + ^Cbh, (33) 

where 



X j ^(d - 2)2(d + 1)252 + 1 [8 + 4(c^ _ a){d - l)d] C7 - ^ [4 - (d - 5)(d - l)d] cg 



(34) 



It should be noted, however, that the validity of Eq. (|34|) has its limitations. Indeed, for small black 
holes the effect of the back reaction should be taken into account and the perturbative approach 
fails to be accurate. Similarly, since the black hole is treated as a system to which a thermal 
description applies, one must require that the change of the Hawking temperature caused by the 
emission of a single quantum of radiation is small. As discussed in Ref [sgl . I40I 1 the condition for 
the thermal description to be self-consistent is 

dT 
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Now, making use of one has 



dT 
dM 



4Gd 
(kod 

a?(hjd 



(2C) 



-d/d-l 



+ 



adujd 



(2C)-('^+2)/(d-i)d(d-2) 



(d - 2)2(4^2 -7d- 4) (2C) 



-(d+4)/(d-l) 



+ 



{d - 4) [4 (d^ - + 2d + 2)c7- {d^ - + 5d-4) dg] (2C) 



-{d+4)/(d-l) 



(36) 



and, therefore, one concludes that mini black holes violate the condition (|H5j) . On the other hand, 
for sufficiently massive black hole the first term in the right hand side of the above equation is 
dominant, and, consequently, the thermal description as well as the perturbative approach is legit- 
imate. Such black holes are, however, thermodynamically unstable and their qualitative behaviour 
is similar to Schwarzschild black hole. Inclusion of the cosmological constant, angular momentum 
or electric charge changes the situation dramatically. 

The entropy of the black hole may be calculated using various method. Two techniques, however, 
are especially well suited for calculations of the entropy in the higher derivative theories. One of 



them is the Wa^ 
redefiniton 43, 



d's Noether charge approach 



41 



42l | whereas the second one is based on the field 



4J]. Here we shall follow the approach, in which S is given in terms of the surface 



integral over the event horizon 

S = 



42 



43 



4^: 



AGd Jh dRabcd ^ 



(37) 



where gab denotes the metric in the subspace orthogonal to the event horizon, and, in the case 
on hand, C is the sum of the Lagrangians given by Eqs. © and Q. It should be noted that to 
calculate the entropy to the required order it suffices to retain in the line element the terms which 
are linear in e. 

The calculation of S consists of three steps. First, it is necessary to express the line element in 
terms of r+, which can be easily achieved by inverting relation (|28() . Simple calculations yield 



C 



lrr + ^(<i-2)(d-l)r 



d-3 

+ 



-{d-2){d- - 1) (6163 + 6263 + 362) ^ 



d-5 



1 



1 

4^ 



{d-l)c^ [2 + d{d - 1){A - hd)] r 



d-5 



+ (d- l)c8 [10 - d(l3 - I2d + 5^2)1 

16a 



(38) 



Equally well one may calculate M{r) using the boundary conditions of the first kind H17j) and take 
the limit as r ^ 00. Further, one has to calculate 

dC 



J- 



abed 



dRabed 



(39) 
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remembering that J°-^'^'^ shares all symmetries of the Riemann tensor. Finally, performing simple 
integration (which actually reduces to the multiplication of the result by the factor uJdf"i^): after 
some algebra, one gets the desired result: 

A 

S = + ATTLUdrie {-2hiR + 262V V + ^^^^tr") \h + ^^'^rf^'+e^ {csRiem^ 



+ C6 



(40) 



where ellipsis denote the terms which are omitted as they will not contribute to the result. All 
quantities are to be calculated at the event horizon and summation is assumed over repeated 
indices: p,q = 0,1 (t, r) and i,j = 2...d + 1. 

Now we are ready to calculate the entropy of the black hole described by Eqs. H13I14|1 with 
(|21ll23p . Substituting formulas collected in Appendix into Eq. (|4Ujl one has 



(41) 



, d{d-l) 



d{d + 1)C3 + i(d2 - 1)C6 + M{d - 1)C7 -^{d- ifcs 



As expected, the higher order corrections to the action modify the standard relation between S 
and the area of the event horizon which is valid only for the Einstein gravity. 
The entropy can also be calculated employing the first law of thermodynamics 

M = TdS + ^fiidQ^, (42) 

i 

where fii are the chemical potentials corresponding to the conserved charges Qj. Making use of Eq. 
(gH) one has 



S = J T-^dM + So 



where the integration constant Sq does not depend on r+, but possibly depends on the coupling 
constants and the spacetime dimension. It should be noted that in the present approach it is 
necessary to retain in the line element all the terms proportional to also. After some algebra 
one obtains the expression describing the entropy, which for 5*0 = coincides with the one given 
by Eq. (jni. 

The integration constant can be determined from the physical requirement that the entropy 
vanishes when the horizon radius shrinks to zero [s^, -^o'^ Lovelock theory it has been 



12 



shown that this condition leads to the results which are identical with those obtained within the 
framework of the Euclidean approach. For d > 4 this procedure gives 5*0 = 0. 

The entropy as given by Eq. (|4H) is expressed in terms of the exact location of the event hori- 
zon, r+, and therefore it depends on the full set of the coupling constants. However, according to 
our previous discussion, one can easily reduce their number by a suitable choice of representation. 
Indeed, expressing the entropy in terms of the total mass as seen by a distant observer (or, equiv- 
alently, C) reduces the number of remaining coupling constants to three. To demonstrate this, let 
us substitute r_|_ given by Ea. H28() into Ea. ()41() and retain the terms up to second order in e. After 
some rearrangement, one obtains 

S = ^(2C)'^/('^-i) + -£!^(2C)('^-2)/('^~i)d263 + i!^(2C)('^-4)/('^-i) \cj (2 + 2d- 3d' + d') 



-^C8 (4 + 5d - ed^ + d^) - hjid - 2)2(1 + 2d) 



(44) 



The entropy described by the above equation depends only on 63, C7 and cg as expected. 

Now, let us compare (jUj) with the analogous result constructed by Lu and Wise [2^. Putting 
in (jUI) d = 2 and Gd = G, one has 

27r2 

S = AttGM' + 647r2e63 + (4c4 + cg) e\ (45) 

Inspection of H45|) shows that it contains the term proportional to 63 , which is absent in the Lu-Wise 
paper. This can be easily understood as Lu and Wise ignored the Gauss-Bonnet term, which, in 
four dimensions, does not affect black hole solution of the field equations. It affects, of course, the 
action itself and consequently the entropy, leading to appearance of a constant term in S that is 
independent of ^A. 

Finally, let us restrict values of the coefficients to its Lovelock combinations. After some ma- 
nipulations it could be shown, that (|41|1 reduces to a simple expression 



ai + ^d{d - 1) + ^d{d - 3)(d - 2){d - 1) 



AGnai 

which is identical with the exact result obtained from a general formula 

d ^ 1 

„ _ i^dr^ sr^ nd 



(46) 



3 



n=l 



for m = 3. 
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V. FINAL REMARKS 

So far we have considered the boundary conditions of the second type only. Now, we shall briefly 
examine appropriate solution constructed with the aid of the conditions (|17j) and ip (oo) = 0. Since 
both solutions are not independent one can treat the solution of the one type as the useful check 
of the other. Before we proceed further let us analyse some general features of the function M(r). 
Nature of the problem and our previous analysis suggests that the solution has the form 



M(r) =M{r) + Ci, 



(48) 



with M (oo) = and M (r+) = r'^ ^/2 — Ci, and, consequently, the total mass of the system as 
seen by a distant observer, Ai, can be obtained from 

SttGd 



M (oo) = Ci 



duJd 



-M, 



(49) 



where Ci is expressed, by construction, in terms of the exact location of the event horizon. One 
expects, therefore, that Ci = C, where C is given by Eq. and this equality can be treated 

as a consistency check. Repeating calculations order by order with the new boundary conditions, 
after some algebra it can be shown that the function M (r) can be written as 

(d-2) (d-^)-^. 

2a 



M{r) 



„3d-3 



+ A2 (r) 



,2d+2 



+ Ci, 



where Ci is given by Eq. (|55|) . 



Ai (r) 



{d-i){d + iy 



-2dc3 + ^ (2 - 3d) C6 - 3 (d - 1) C7 + ^ (d - 1) cg 



+ 



id-2){d-l)id+lY 



(36263 + 46163 + 86i) 



(50) 



(51) 



and 



A2 (r) 



(d-1) 



C3 



d{d+l) {3d + 5) - ^ (d + 1) (3 - 3(i - 2d^) 



a 

-— (10 + 16d - 21d^ - 7d^) + — (2 + 25d - 24d2 
1 



7d^ 



^ (d - 2) (d - 1) [26263 (d + 1) (d + 2) + 6163 (d + 1) (3d + 5) 

+bl (9 + 19d + 4d2)] . 



(52) 
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Similar calculations carried out for the function ^l^{r) yield 
2 



2d-2 



(d - 2) (d - 1) (d + 1) [6263 {d + 2) + bibs {2d + 3) + 6i {2d + 5) 



p2d-\-2 



{d-i){d + i) 



C3d (3 + 2d) - C6 (1 - 2(i - d^) - -C7 {2 -2d- d^) 



--C8 {2 -3d- d^) 



r2d+2 ■ 



(53) 



Since, by assumption, the radius of the event horizon is treated as the exact quantity now, the 
thus derived line element may be easily employed in construction of the entropy. First, observe 
that the Hawking temperature calculated with the aid of the Eq. (|3U() is given in terms of and 
depends on all relevant coefficients. On the other hand the relation H31|) is independent of bi, 62, C4 
and C5. This behaviour can be ascribed to the possibility to remove the dependence of the line 
element on this very coefficients at the expense of modifications of the equations of motion of test 
particles. Indeed, it could be demonstrated that by means of the field redefinition of the form 



(54) 



where 



Ail^ =q\"R9ab + q'2'Rab 



(1) 



(55) 



and 



^ab — Ql ^ 9ab + Q2 ^^ab + % 9ab^cdef^ + <?4 9ab^cd^ 
+ (li'' RacRl + RacdeRb'^'^'^ ^ 



ycd 



(56) 



one can remove all the terms in the action except these proportional to the parameters 63, C7 and 
cg. The coefficients and q^^^ can be determined by solving, at each order of the expansion, 
the appropriate systems of equations. As the result of the field redefinition 1)541) . one obtains two 
additional terms ROR and Rab^R"''', which can also be removed from the action functional. It 
should be noted, however, that such terms appear naturally in the effective action of the quantized 



massive fields in a large mass limit 



45 



4S 



The geodesic equation after the field redefinition becomes 



d^x^ 



+ ds ds 



(1) 

mk;j 



1.(1) \dx^(kf_ 
2 i^-^'^j ds ds 



2 im 



e 9 



,(2) _1,(2) \d^d^_ 2,{l)i 

^mk;j 2 i'^;'"/ ds ds 



A 



(1) 

mk;j 



2 jk;m 



dx^ dx^ 
ds ds 



0, (57) 
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where g is the affine parameter of the original metric. Now, one can repeat the arguments of 
Ref. 20]. Both A4 and T can be measured at infinity and do not depend on the particular 
form of the equations of motion. Consequently, the temperature mass relation is independent of 
the removed terms. However, to determine the radius of the event horizon one performs local 
measurements and the equations of motion of test particles are important. 

It could be easily seen that, as expected, the entropy is precisely the same for both choices of 
boundary conditions and is described by the formula H41|). Since the calculations for both types 
of the boundary conditions have been carried out independently, this equality may be regarded as 
the additional consistency check. 

APPENDIX 

In this Appendix we list some formulas used in Section IIVI The expansions of the components 
of the Riemann tensor and some of its contractions up to the first order in e, which are necessary 
in calculation of the entropy, are given by 



= -T - '5i«J ) + O ■ (A.2) 



and 



, , <'ii_21 _ - 2) 1) d(d + 1) ^ ^ ^^ jj 

2rl 2arl 



j,r^_^,,^(A^m±:J^^M±Il + o{e-) (A.4) 

2ari 



2arl 



where k,l = 2, d + 1. 
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